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Abstract 



In this paper we are interested in establishing stability estimates in the 
inverse problem of determining on a compact Riemannian manifold the elec- 
tric potential or the conformal factor in a Schrodinger equation with Dirich- 
let data from measured Neumann boundary observations. This information 
is enclosed in the dynamical Dirichlet-to-Neumann map associated to the 
Schrodinger equation. We prove in dimension n > 2 that the knowledge of 
the Dirichlet-to-Neumann map for the Schrodinger equation uniquely deter- 
mines the electric potential and we establish Holder-type stability estimates 
in determining the potential. We prove similar results for the determination 
of a conformal factor close to 1. 

Keywords: Stability estimates, Schrodinger inverse problem, Dirichlet-to- 
Neumann map. 
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1 Introduction and main results 



This paper is devoted to the study of the following inverse boundary value prob- 
lem: given a Riemannian manifold with boundary determine the potential or the 
conformal factor of the metric in a dynamical Schrodinger equation from the ob- 
servations made at the boundary. Let {Ai , g) be a compact Riemannian manifold 
with boundary dAi. All manifolds will be assumed smooth (which means C°°) 
and oriented. We denote by Ag the Laplace-Beltrami operator associated to the 
metric g. In local coordinates, g(x) = {gjk), Ag is given by 




(1.1) 



Here (g-''^) is the inverse of the metric g and detg = det(gjfc). Let us consider 
the following initial boundary value problem for the Schrodinger equation with 
bounded electric potential q G L°°{M-) 

' {idt + + q{x)) u = 0, in {0,T) x M 
< ti(0,-) =0, in M (1.2) 

^ u = f, on (0, T) X dM 
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where / G H^{{0, T) x dJ\4). Denote by v the outward normal vector field along 

n 

the boundary dM., so that g^^VjUk = 1. Further, we may define the dynami- 

j,k=i 

cal Dirichlet-to-Neumann map Ag^ g associated to the Schrodinger equation by 



du 

j,k=i 



(1.3) 

(0,T)x9A4 



Unique determination of the metric g = (gjfc) from the knowledge of the Dirichlet- 
to-Neumann map Ag g is hopeless: as was noted in [33 J in the case of the wave 
equation, the Dirichlet-to-Neumann map is invariant under a gauge transformation 
of the metric g. Namely, if one pulls back the metric g by a diffeomorphism ^ : 
— > which is the identity on the boundary "^\dM = W into a new metric 
^*g, one has A,j,«g ^ = Ag g. The inverse problem has therefore to be formulated 
modulo the gauge invariance. However, we will restrict our inverse problem to a 
conformal class of metrics (for which there is no gauge invariance): knowing A^g^g, 
can one determine the conformal factor c and the potential q? 

In the case of the Schrodinger equation, Avdonin and Belishev gave an affirma- 
tive answer to this question for smooth metrics conformal to the Euclidean metric 
in ||3l. Their approach is based on the boundary control method introduced by Beli- 
shev im and uses in an essential way a unique continuation property. Because of the 
use of this qualitative property, it seems unlikely that the boundary control method 
would provide accurate stability estimates. More precisely, when 7W is a bounded 
domain of M", and g, q £ C^(A^) are real functions, Avdonin and Belishev lO 
show that for any fixed T > the response operator (or the Neumann-to-Dirichlet 
map) of the Schrodinger equation {igdtu + Au — qu) = uniquely determines the 
coefficients g and q. The problem is reduced to recovering g, q from the boundary 
spectral data. The spectral data are extracted from the response operator by the use 
of a variational principle. 

The uniqueness in the determination of a time-dependent electromagnetic po- 
tential, appearing in a Schrodinger equation in a domain with obstacles, from the 
Dirichlet-to-Neumann map was proved by Eskin fTV'l. The main ingredient in his 
proof is the construction of geometrical optics solutions. In [2], Avdonin, Lenhart 
and Protopopescu use the so-called BC (boundary control) method to prove that the 
Dirichlet-to-Neumann map determines the time-independent electrical potential in 
a one dimensional Schrodinger equation. 

The analogue problem for the wave equation has a long history. Unique deter- 
mination of the metric goes back to Belishev and Kurylev 1 6 1 using the boundary 
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control method and involves works of Katchlov, Kurylev and Lassas ll25l . Kurylev 
and Lassas lITTl ) and Anderson, Katchalov, Kurylev, Lassas and Taylor HI. In 
fact, Katchalov, Kurylev, Lassas and Mandache proved that the determination of 
the metric from the Dirichlet-to-Neumann map was equivalent for the wave and 
Schrodinger equations (as well as other related inverse problems) in 1261 . Identi- 
afiability of the potential was proved by Rakesh and Symes 1301 in the Euclidian 
case (g = e) using complex geometrical optics solutions concentrating near lines 
with any direction w G to prove that Ag^g determines q{x) uniquely in the 

wave equation. This result was generalized by Ramm and Sjostrand [31] and Eskin 
lITSl [T9l to the case of q depending on space and time. Isakov 1221 also considered 
the simultaneous determination of a potential and a damping coefficient. 

As for the stability of the wave equation in the Euclidian case, we also refer 
to |[35l and 1231; in those papers, the Dirichlet-to-Neumann map was considered 
on the whole boundary. Isakov and Sun [23] proved that the difference in some 
subdomain of two coefficients is estimated by an operator norm of the difference 
of the corresponding local Dirichlet-to-Neumann maps, and that the estimate is 
of Holder type. Bellassoued, Jellali and Yamamoto 1 10] considered the inverse 
problem of recovering a time independent potential in the hyperbolic equation from 
the partial Dirichlet-to-Neumann map. They proved a logarithm stability estimate. 
Moreover in [29] it is proved that if an unknown coefficient belongs to a given 
finite dimensional vector space, then the uniqueness follows by a finite number 
of measurements on the whole boundary. In [[Tl, Bellassoued and Benjoud used 
complex geometrical optics solutions concentring near lines in any direction to 
prove that the Dirichlet-to-Neumann map determines uniquely the magnetic field 
induced by a magnetic potential in a magnetic wave equation. 

In the case of the anisotropic wave equation, the problem of establishing sta- 
bility estimates in determining the metric was studied by Stefanov and Uhlmann in 
ll33l l34l for metrics close to Euclidean and generic simple metrics. In a previous 
paper (TT), the authors also proved stability estimates for the wave equation in de- 
termining a conformal factor close to 1 and time independent potentials in simple 
geometries. We refer to this paper for a longer bibliography in the case of the wave 
equation. 

The inverse problem for the (dynamical) Schrodinger equation seems to have 
been a little bit less studied. In the Euclidean case, there are extensive results by 
Bellassoued and Choulli |8| where a Lipschitz stability estimate was proven for 
time independent magnetic potentials. The stability problem in determining a time 
independent potential in a Schrodinger equation from a single boundary measure- 
ment was studied by Baudouin and Puel [4 j. They established Lipschitz stability 
estimates by a method based essentially on an appropriate Carleman inequality. In 
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the above mentionned papers, the main assumption is that the part of the boundary 
where the measurement is made must satisfy a geometrical condition (related to ge- 
ometric optics condition insuring observality). Recently, Bellassoued and ChouUi 
showed in tOJ that this geometric condition can be relaxed provided that the poten- 
tial is known near the boundary. The key idea was the following : the authors used 
an FBI transform to change the Schrodinger equation near the boundary into a heat 
equation for which one can use a useful Carleman inequality involving a boundary 
term and without any geometric condition. 

The main goal of this paper is to study the stability of the inverse problem for 
the dynamical anisotropic Schrodinger equation. We follow the same strategy as in 
ifTTI inspired by the works of Dos Santos Ferreira, Kenig, Salo and Uhlmann 1161 . 
Stefanov and Uhlmann ll33l[34l and Bellassoued and ChouUi HI. 



1.1 Weak solutions of the Schrodinger equation 

First, we will consider the initial-boundary value problem for the Schrodinger 
equation on a manifold with boundary ( 11.21 ). This initial boundary value prob- 
lem corresponds to an elliptic operator — Ag given by (11.11 ). We will develop an 
invariant approach to prove existence and uniqueness of solutions and to study their 
regularity proprieties. 

Before stating our first main result, we recall the following preliminaries. We 
refer to 1 24 1 for the differential calculus of tensor fields on a Riemannian manifold. 
Let {A4 , g) be an n-dimensional, n > 2, compact Riemannian manifold, with 
smooth boundary and smooth metric g. Fix a coordinate system x = [xi, . . . , Xn] 



and let 



a a 

axi ' • ■ ■ ' dXn 



be the corresponding tangent vector fields. For x £ Ai, the 



inner product and the norm on the tangent space Tx^A are given by 

n 

g{X,Y) = {X,Y)^= J2 SjkajPk, 
j,k=i 

i=i * 1=1 * 

If / is a function on A4, the gradient of / is the vector field Vg/ such that 

X(/) = (Vg/,X)g 

for all vector fields X on A4. This reads in coordinates 

t,j=i ■' 
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1 /2 

The metric tensor g induces the Riemannian volume dvg = (det g) ' dxi A • • • A 
dxn- We denote by L'^{M.) the completion of C°°{M) endowed with the usual 
inner product 

(/i,/2)=/ /i(x)^Mdv-, h,f2^C^{M). 
Jm 

The Sobolev space H^{Ai) is the completion of C°^(7W) with respect to the norm 

II ■ ll_ffi(M)' 

WfWn^M) = II/IIl2(>!) + ll^g/llL2(A4) • 
The normal derivative is 

d.u:=V,u.i.= Eg'S^ (1-5) 

j,k=l 

where u is the unit outward vector field to dAi. Moreover, using covariant deriva- 
tives (see 1.20,1 ). it is possible to define coordinate invariant norms in H^{M.), 
k>0. 

Before stating our main results on the inverse problem, our interest will focus 
on the study of the initial boundary problem (II. 2I ). when u is a weak solution in 
the class C{0,T; H^{M)) D C^{0,T; H-^{M)). The following theorem gives 
conditions on / and q, which guarantee uniqueness and continuous dependence on 
the data of the solutions of the Schrodinger equation (11.21 ) with non-homogenous 
Dirichlet boundary condition. 

Theorem I LetT > be given. Suppose that f G H'^{{0,T) x dM) and q £ 
W^'°°{A4). Then the unique solution u of il.2\) satisfies 

u G C{0,T;H^{M))nC^{0,T-H'^{M)), (1.6) 

duue L'^{{0,T) xdM). (1.7) 

Furthermore, there is a constant C = C{T, Ai) > such that 

ll^!^"llL2((o,T)xaA^) ^ ^ \\f\\m{{o,T)xdM) ■ (1-^) 

The Dirichlet-to-Neumann map Ag^g defined by (fJ.3l) is therefore continuous and 
we denote by \\Ag,q\\ its norm in c(h^{{0, T) x dM),L'^{{0,T) x dM)). 

Theorem [T] gives a rather comprehensive treatment of the regularity problem for 
(11.21 ) with stronger boundary condition /. Moreover, our treatment clearly shows 
that a regularity for / G H^{{0,T) x dM.) is sufficient to obtain the desired 
interior regularity of u on (0, T) x M while the full strength of the assumption 
/ G i/^((0,T) X dM) is used to obtain the desired boundary regularity for dyU 
and then the continuity of the Dirichlet-to-Neumann map Ag g. 
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1.2 Stable determination 



In this section we state the main stability results. Let us first introduce the admis- 
sible class of manifolds for which we can prove uniqueness and stability results in 
our inverse problem. For this we need the notion of simple manifolds 1341 . 

Let {A4 , g) be a Riemannian manifold with boundary ^^A , we denote by D 
the Levi-Civita connection on {M-, g). For a point x G dA4, the second quadratic 
form of the boundary 

is defined on the space Tx{dM). We say that the boundary is strictly convex if the 
form is positive-definite for all x G dM.. 

Definition 1 We say that the Riemannian manifold {M , g) (or that the metric g) 
is simple in Ai, if dAi is strictly convex with respect to g, and for any x £ Ai, 
the exponential map exp^ : exp~^(7W) — > Ai is a diffeomorphism. The latter 
means that every two points x;y £ Ai are joined by a unique geodesic smoothly 
depending on x and y. 

Note that if {Ai , g) is simple, one can extend it to a simple manifold Ai i such that 
Ml DM. 

Let us now introduce the admissible set of potentials q and the admissible set 
of conformal factors c. Let Mq > 0, > 1 and e > be given, set 

^(Mo) = [qe W'^'^iM), |kllvKi-{M) < ^^o} , (1-9) 

and 

^(Mo,A;,e) = 

{c G C°°{M), c > inM, 111 - c||ci(^) < e, \\c\\ck(^M) ^ ^o} • (1-10) 
The main results of this paper are as follows. 

Tlieorem 2 Let [M , g) be a simple compact Riemannian manifold with boundary 
of dimension n > 2 and let T > 0. There exist constants C > and s G (0, 1) 
such that for any gi, ^2 S ■S{Mq), qi = q2 on dM, we have 

-^2||l2(^) < C||Ag,g, - Ag^g^f (1.11) 

where C depends on M, T, Mq, n, a and s. 
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By Theorem [21 we can readily derive the following uniqueness result 

Corollary 1 Assume that T > 0. Let qi, q2 G ^(Mq), qi = q2 on dM. Then 
Ag^qj = Ag implies qi = q2 everywhere in M.. 

Theorem 3 Let {Ai , g) be a simple compact Riemannian manifold with boundary 
of dimension n > 2 and let T > 0. There exist k>l, e>0, 0<s<l and 
C > such that for any c G "^(Mq, k, e) with c = 1 near the boundary dM., the 
following estimate holds true 

||l-c||^.(^) <C||Ag-Aeg|r (1.12) 

where C depends on M., Mq, n, e, k and s. 

By Theorem [3l we can readily derive the following uniqueness result 

Corollary 2 Let (Ai ,g) be a simple compact Riemannian manifold with boundary 
of dimension n > 2 and let T > 0. There exist k > 1, e > 0, such that for any 
c G 'Tf{Mo,k,e) with c = 1 near the boundary dM, we have that Acg = Ag 
implies c = 1 everywhere in M. 

Our proof is inspired by techniques used by Stefanov and Uhlmann fi^, and 
Dos Santos Ferreira-Kenig-Salo-Uhlmann 1 16] which prove uniqueness theorems 
for an inverse problem related to an elliptic equation. Their idea in turn goes back 
to the pioneering work of Calderon |[T3l . We also refer to Bukhgeim and Uhlamnn 
|[r2il . Cheng and Yamamoto fTSl, Hech-Wang ||2TI and Uhlmann B6l as a survey. 

The outline of the paper is as follows. In section 2 we collect some of the 
formulas needed in the paper. In section 3 we study the Cauchy problem for the 
Schrodinger equation and we prove Theorem [U In section 4 we construct special 
geometrical optic solutions to Schrodinger equations. In section 5 and 6, we estab- 
lish stability estimates for related integrals over geodesies crossing M. and prove 
our main results. 



2 Geodesical ray transform on a simple manifold 

In this section we first collect some formulas needed in the rest of this paper and 
introduce the geodesical ray transform. Denote by divX the divergence of a vector 
field X G H^{TM) on M, i.e. in local coordinates, 

^ n ^ d 
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If X G H^(TAi) the divergence formula reads 

/ divXdv"=/ (X,i/)da^-^ (2.2) 
JM JdM 

and for / G H^{Ai) Green's formula reads 

/ divX/dv^ = -/" {X,VJ) dv2+ [ (X,z.)/d^7^-i. (2.3) 
Jm Jm ^ JdM 

Then if / G H^{M) and w G H'^{M), the following identity holds 

/ A^wfdv^ = - [ {V^w,VJ) dv^+ [ d.wfda^-'. (2.4) 
Jm Jm Jsm 

Let u G {M) and be a smooth real vector field. The following identity holds 
true (see |[37l ) 

<Vg^,Vg((iV,VgtJ)g)>^ 

= DN{V^v,V^v) + ^div(lVg^l^iV) - ^ |Vg^;|^div(iV) (2.5) 

where D is the Levi-Civita connection. 

For X ^ M and 9 G Ta;7W we denote by 7^. g the unique geodesic starting at the 
point X in the direction 9. We consider 

SM = \{x,9) eTM; \9\^ = l} , 

S*M = {ix,p)eT*M; Ng = l} 

the sphere bundle and co-sphere bundle of A4. The exponential map exp^ : 
TxM — > M is given by 

^^Vx{v) = lx,v{\'^\gv) = -ix,v{rv), r=\v\^. (2.6) 

A compact Riemannian manifold (TW , g) with boundary is called a convex non- 
trapping manifold, if it satisfies two conditions: 

(a) the boundary dM^ is strictly convex, i.e., the second fundamental form of the 
boundary is positive definite at every boundary point, 

(b) for every point x ^ M and every vector 9 G T^M., 0/0, the maximal 
geodesic ^x,e{t) satisfying the initial conditions 7x,e(0) = x and 7x,e(0) = 9 
is defined on a finite segment [r_(x, 9), t+(x, 9)]. We recall that a geodesic 
7 : [a,b] — > M is maximal if it cannot be extended to a segment [a — ei , 6+ 
£2], where > and ei + 62 > 0. 
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The second condition is equivalent to all geodesies having finite length in M. . 

An important subclass of convex non-trapping manifolds are simple manifolds. 
We way that a compact Riemannian manifold {Ai , g) is simple if it satisfies the 
following properties 

(a) the boundary is strictly convex, 

(b) there are no conjugate points on any geodesic. 

A simple n-dimensional Riemannian manifold is diffeomorphic to a closed ball in 
M", and any pair of points in the manifold are joined by an unique geodesic. 

Now, we introduce the submanifolds of inner and outer vectors of SA4 

d±SM = {{x,e) G SM, X G dM, ±{e,v{x)) < 0} (2.7) 

where u is the unit outer normal to the boundary. Note that d-^-SAd and d-SA4 
are compact manifolds with the same boundary S{dA4), and dSAi = d+SJv[ U 
d^SM. For {x,e) E d+SM, we denote by 73,^0 : [0,r+(x,6')] — > M the 
maximal geodesic satisfying the initial conditions 7a:,6»(0) = x and '^x,e{^) = ^■ 
Let C°°{d+SM.) be the space of smooth functions on the manifold d+SM.. The 
ray transform (also called geodesic X-ray transform) on a convex non trapping 
manifold M. is the linear operator 

Z : C°°{M) C°°{d+SM) (2.8) 

defined by the equality 

X/(x, 9) = r'^^^'^ f{lx,e{t)) dt. (2.9) 
Jo 

The right-hand side of ( 12.91 ) is a smooth function on dj^SJvi because the integration 
limit T+(x, 6) is a smooth function on d+SM, see Lemma 4.1.1 of |f32l. The ray 
transform on a convex non trapping manifold M can be extended as a bounded 
operator 

X : H^{M) H^{d+SM) (2.10) 
for every integer k>l, see Theorem 4.2.1 of 021 . 

The Riemannian scalar product on T^M induces the volume form on Sx-M., 
denoted by dujx{9) and given by 

n 

dw^(0) = ^{-ife^de"" ^de^ ^ dr. 

k=l 
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We introduce the volume form dVg" on the manifold SM by 

dv2"-i(x,0) = |d^^,(0) A dv^l 

where dVg is the Riemannnian volume form on M.. By Liouville's theorem, the 
form dVg""^ is preserved by the geodesic flow. The corresponding volume form 
on the boundary dSM = {(x, 6) G SM^ x G dM} is given by 

where dUg"^ is the volume form of ^^A. 

Let L'^^{d+SM.) be the space of square integrable functions with respect to 
the measure ^{x, 0) d(7g"~^ with ^{x, 6) = \ {6, This real Hilbert space is 

endowed with the scalar product 

{u,v)L2^a+SM)= f u{x,e)vix,9)f,{x,e)dal''-^. (2.11) 
^ Jd+SM 

The ray transform X is a bounded operator from L^(A^) into L'^{^-^-SAi). The 
adjoint I* : Ll{d+SM) L^{M) is given by 

x*i){x)= f ^*{x,e)duj^{e) (2.12) 

where ip* i?, the extension of the function ip from d^SM to SM constant on every 
orbit of the geodesic flow, i.e. 

V'*(x,0)=V(7x,e(r+(x,0))). 

Let {M. , g) be a simple metric, we assume that g extends smoothly as a simple 
metric on A^i ^ M. Then there exist Ci > 0, C2 > such that 

Ci < ||X*X(/)||^,(^^) < C2 (2.13) 

for any / E L^(A^).If V is an open set of the simple Riemannian manifold 
{M. 1 , g) , the normal operator X*X is an elliptic pseudodifferential operator of order 
— 1 on y whose principal symbol is a multiple of |^|g (see ll34l ). Therefore there 
exists a constant Ck > such that for all / G H^{V) compactly supported in V 

U*Af)\\m+^iM,)<Ck\\f\\miyy (2.14) 
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3 The Cauchy problem for the Schrodinger equation 



In this section we will establish existence, uniqueness and continuous dependence 
on the data of solutions to the Schrodinger equation (11.21 ) with non-homogenous 
Dirichlet boundary condition / G H^{{0, T) x dM). We will use the method of 
transposition, or adjoint isomorphism of equations, and we shall solve the case of 
non-homogenous Dirichlet boundary conditions under stronger assumptions on the 
data than those in [4]. 

Let us first review the classical well-posedness results for the Schrodinger 
equation with homogenous boundary conditions. After applying the transposition 
method, we establish Theorem [T] 

3.1 Homogenous boundary condition 

Let us consider the following initial and homogenous boundary value problem for 
the Schrodinger equation 

' {idt + + q{x))v{t,x) = F{t,x) in (0,r)xX, 
v{<d,x) =0 in A^, (3.1) 

v{t,x) =0 on (0,r) X dM. 

Firstly, it is well known that if F G -L^(0, T; Lp'{M)) then ( 13.11 ) admits an unique 
weak solution 

V ^C{Q,T-L^{M)) . (3.2) 

If we multiply both sides of the first equation ( 13.11 ) by v and integrate over M., we 
obtain 



Im 



idtv{t)vdYl- / \V^v{t,x)\l+q\v{t,x)\'' dwl 
M Jm ^ . 



1 d 
2dt 



dv^ = Im / (F{t,x)v{t,x)+q\v{t,x)\'^)dv'^. (3.3) 
takeao(t) = lb(i)llL2(^)> * ^ (0,T), weget 

^(ag(t)) <c(||F(t,.)IL2(^)«o(t) + a2(t)), Vte(0,r), (3.4) 
which implies that a'^it) < C {\\F{t, + cto(i)) ^^id 

\\v{t)\\L,^j^^<CT l|i^(i,-)llL2(^) dt, VtG(0,r). (3.5) 
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Lemma 3.1 LetT >Oandq£ Suppose that F e L^{0,T; H^{M)). 

Then the unique solution v of ( I3.1I ) satisfies 



veC{0,T;H^{M)). (3.6) 

Furthermore there is a constant C > such that for any F G -Z>^(0, T; Hq{M)), 
we have 

')\\hI{M) - ^ II-^IIli(0,T;H1(A1)) ■ (^•'^) 

Proof . Using the classical result of existence and uniqueness of weak solutions in 
Cazenave and Haraux lfT4l (set for abstract evolution equations), we obtain 

veC{0,T;H^{M)). (3.8) 

Multiplying the first equation of (13.11 ) by AgU and using Green's formula, we get 

Im / idtv{t)Agvdv2 - \Agv{t)\^ + qvAgV dv\ 
_Jm Jm 

= Im^(VgF(t,x),Vg^) dv«. 

Letai(t) = ||Vgv(t)||^2(^),t G (0,T). Then, by (ESll, we have 
d 



(3.9) 



^^(a?(t)) <C(||F(t,.)||^^i(A4)«iW+«?Wj, ViG(0,r), (3.10) 

which implies that a']^(t) < C {\\F{t,-)\\^i^j^-^ +ai(t)^ and by Gronwall's lemma 
we find 

\HmHl[M)<CT I dt, VtG(0,r). (3.11) 



The proof of (13.71 ) is complete. □ 



Lemma 3.2 LetT >0 and q G L°°{M). Suppose that F G W^'^{0,T; L'^{M)) 
such that F{0, •) = 0. Then the unique solution v of (13.11 ) satisfies 

V G C\0,T; L\M)) nC{0,T; H^{M) n H^{M)). (3.12) 

Furthermore there is a constant C > such that for any 7] > small, we have 
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Proof . If we consider the equation satisfied by dtv, (13.21 1 provides tiie following 
regularity 

V £C^{id,T-L^{M)). 

Furthermore by ( I3.5I ). there is a constant C > such that the following estimate 
holds true 

\\dtv{t,-)\\L,^M^ < c£\\dtF{t,.)\\^,^^^ dt yt e (0,T). (3.14) 

Then, by dlJ]), we see that A^v = -idtv + F £ C{0,T; L^{M)) and therefore 

V G C{0,T;H\M)). 
This complete the proof of ( 13.121 ). 

Next, multiplying the first equation of (13.11 ) by v and integrating by parts, we obtain 



Re 



idtvit,x)vit,x)dy^- / (\Vgv{t)\^ -q\v\^) dv^ 



Re/ F(t,x)u(t,x)dv" 
JM 

Re y [J^ dtF{s, x) ds^ v{t, x) dv^. (3. 15) 



Take now ai{t) = ||Vgt;(t)||^2(^)- Then there exists a constant C > such that 
the following estimate holds true 



ajit) < C 



+ 



\\dtv{tr)\\L'2(M) II^(^'-)IIl2(A4) + II^(0IIl2(A1) 



\v{t,x)dtF{s,x)\ dv^ds 



'0 JM 

Using (13.141 ) and ( |331 ). we get 



(3.16) 



Q?(t) < C 



l^*-^llLi{0,r;L2(A4)) II-^IIl1(0,T;L2(A4)) 



Li(0,T;L2(7K)) 



Thus, we deduce (13.131 ). and this concludes the proof of Lemma 1X2] 



(3.17) 

□ 



Lemma 3.3 LetT>0,qe W^'°°{M) be given and let M' = Li(0, T; Hl{M)) 
or .j^ = (0, T; Lp'{M.)). Then the mapping F i— )• dyV where v is the unique so- 
lution to (13.11 ) is linear and continuous from to L'^{{0, T) x dM). Furthermore, 
there is a constant C > such that 



l^i''"llL2((o,T)x9A4) - ^ ll-^lljf 



(3.18) 
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Proof . Let be a vector field on M. such that 



N{x) = u{x), x£dM; |A^(x)| < 1, x £ M. 



(3.19) 



Multiply both sides of the first equation in (13.11 ) by {N, Vgw)^ and integrate over 
(0,r) X X, this gives 



JM 



F{t,x) (A^,VgU)g dv^dt 



idtv {N, VgU)g dv^ dt + / /" Ag {N, VgU)g dv^ dt 

Jo J M 

q{x)v {N ,V £v) ^ dv^dt 



/O JM 



/O JM 

h + h + h 



(3.20) 



Consider the first term on left side of (13.20b : integrating by parts with respect t, we 
get 



JM 



M 



idtv {N, VgU)g dv^ dt 
1 ^ 

^^(^,VgI;)g dv^ -i 



(3.21) 



JM 



v{N,Vgdtv)^ dv^dt 



T 



i I viT, x) {N, V^v{T, dv^ - i / / {N, Vg{vdtv))^ dv^ dt 



M 



JM 



T 



+ 1 



dtlj{N,Vgv)^ dv"dt. 



'0 JM 

Then, by (12.3b . we obtain 

f-T 

idtv{N,Vgv)^dv2dt 



Re 



JM 



i [ v{T,x){N,Vgv{T,x)) dv'l + i [ [ div{N)vdtvdv2dt 
JM Jo JM 



— z 



vdtvda'^-'dt 



JdM 



ij viT,x){N,V^viT,x))^dv2+ I I (VgI;,Vg(div(iV)i;))g dv^dt 



M 



JM 



+ 



[ [ Fdiv{N)vdv^dt- [ [ gdiv(iV) dv^ dt 
Jo Jm Jo Jm 
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Jo JdM 



JdM 



The last term vanishes, using (13.131 ) or (13.71 ). we conclude that 

|Re/i|<C||F||2^. 
On the other hand, by Green's theorem, we get 



(3.22) 



l2 = U^A^v{N,V^v)^ dv«dt 



T 



V^v,V,i{N,V,v)^)) dv^dt + 



'0 JM 

Thus by (12.51 ). we deduce 



JdM 



\d,v\' da^-'dt-- 
lO JdM ^ 



/O JdM 



|Vgt;|2 dal'^dt 



+ 

Using the fact 



/ / B^N{V^v,V^v)dvldt-\ \ [ |VgZ;|2div(iV)dVg"dt. 
Jo JM ^ Jo J M 



where V,- is the tangential gradient on dA4 , we get 



Re/2 = - 
2 



'0 JdM 
(■T 



\ff |VgHgdiv(Ar)dv^dt. 

^ Jo JM 

Finally by Lemma [37T] and 1X2] we have 

|Re/3|<||F||2^. 
Collecting (Ill4l ). (Ill3] ). (IHIl ) and (IHOl ). we obtain 



/ / L>giV(Vg?;, Vg?;) dv" dt 



T 



/O 

This completes the proof of (13.181 ). 



(3.23) 
(3.24) 

(3.25) 

□ 
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3.2 Non-homogenous boundary condition 

We now turn to the non-homogenous Schrodinger problem (11.21) . Let 

^ = L\0,T;H^{M)) or = H^{0,T; L"^ {M)). 

By (•, •)^, ^, we denote the dual pairing between J^f" and J^. 

Definition 2 LetT > 0, q e W'^^°°{M) and f € ^^((o, T) x dM), we say that 
u G M" is a solution of (11.21 ) in the transposition sense if for any F G .3^ we have 

(u,F)^, ^= [ [ f{t,x)d^v{t,x)da':-Ut (3.26) 
Jo JdM 

where v = v{t, x) is the solution of the homogenous boundary value problem 
' {idt + Ag + q{x))v{t,x) = F{t,x) in {0,T) x M, 

< v{T,x) = in M, (3.27) 

v{t,x) =0 on (0,r) X dM. 

One gets the following lemma. 

Lemma 3.4 Let f G -Z>^((0, T) x dAi). There exists a unique solution 

u G C(0, T; H-^{M)) n i?"^(0, T; L^{M)) (3.28) 
defined by transposition, of the problem 

' {idt + \ + q{x))u{t,x) =^ in (0,r)x7W, 

< u(0,x) = in M, (3.29) 

u{t,x) = f{t,x) on {0,T)xdM. 

Furthermore, there is a constant C > such that 

ll^llc(0,T;//-i(A1)) + ll^ll//-i(0,T;L2(A4)) - ^ II / II L2((0,r) x9A4) ■ (3.30) 

Proof. Let F G ^ = L^{0,T; H^{M)) or Jif = H^{0,T; L'^{M)). Let 
V G C{0,T; Hq{A4)) solution of the backward boundary value problem for the 
Schrodinger equation (13.271 ). By Lemma [331 the mapping F i— )• |^ is linear and 
continuous from to F^((0, T) x M) and there exists C > such that 

ll^llc(o,T;Hi(A4)) - C'll^ll^ (3.31) 
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and 

\\^'^'"\\L^{{0,T)xdM) - ^ W^W.^ ■ (3.32) 
We define a linear functional I on the linear space M' as follows: 

l{F) = {f,d,v)^ 

where v solves (I3.27I ). By ( I3.32I ). we obtain 

\m\ < 11/11 

It is known that any Unear bounded functional on the space J^f can be written as 

where u is some element from the space J^'. Thus the system (13.291 ) admits a 
solution u G J^' in the transposition sense, which satisfies 



This completes the proof of the Lemma. □ 

In what follows, we will need the following estimate for non-homogenous elliptic 
boundary value problem. 

Lemma 3.5 Let ^ e H~^{M) and 4> G H^{dM). Let w £ H^{M) the unique 
solution of the following boundary value problem 



(3.33) 



/S.^w{x) = ip{x) in M., 

w{x) = (j) on dAi, 

then the following estimate holds true 

(3.34) 

Proof . We decompose the solution w of (13.331 ) as w = wi + W2 with wi and W2, 
respectively, solution of 

l\gWi{x) = ii{x) in M, ( Agty2(x) = in M, 

■ I (3.35) 
wi{x)=Q on dM, \ W2{x) = (j) on dM, 

Since — Ag is an isomorphism from //q (A^) to H^^{Ai), we have 

Ikill^i(^) < C7||V'||^_i(_^) . (3.36) 
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Next, it is well known that (see fZm ) 

\\w2\\l2{M) ^ C 



Now, we shall show that 



< C 



(3.37) 



(3.38) 



Indeed, let h G H^{dM) and 9 € C^(0,r), 6 > Let v solve the following 
initial boundary value problem for the wave equation 

' Ag)v(t,x) = in (0,T)x7W, 

u(0,x) = ajv(0,x) = in A^, (3.39) 

^ v{t,x) = h{x)9{t) on (0,T) x dM. 

Then we have (see 1251) 

V G C{0,T;H\M))nC\0,T;L'^{M)). 
Furthermore there exist C > such that 

ll^llc(0,T;//i(A1)) + ll^llci(0,T;L2(A4)) + W^'^'^W L2{{0,T)xdM) 

<C\\h\\j,,^gj^y (3.40) 



Multiplying both sides of ( 13.391 ) by W2 and integrating over (0, T) x 7W, we get 

= [ [ {dfv - Agv)w2{x)dv':dt 
Jo Jm 



M 



dtv(T,x)w2ix)dv'^- I I duvcpix) da2~^ dt 



JdM 



+ / e{t)dt [ h{x)d^W2da 
Jo JdM 

Then, by (l340l) and (l337]) . one gets 



n-l 

g 



(3.41) 



/ 

JdM 



h{x)d^W2da'2 ^ 



< C'( ||^2((o,T)x9A4) W\\L'^{dM) 
+ II''^2|Il2(^) II^IIc1(0,T;L2(A1)) ) 
Il2{9A4) W^H^dM) ■ 



< c 



(3.42) 
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This implies 
Furthermore, Green's formula yields 

/ |Vgu;2p dv^' < \\duW2\\H-i^QM) UWn^idM) 

J M 

iVxom ( 1X441 ) and ( 13371 ). we get 

Both ( 1X451 ) and ( 13361 ) yield likewise 

M\h-\m) + 

This completes the proof of (I3.34! ). 



< c 



(3.43) 
(3.44) 

(3.45) 

(3.46) 

□ 



3.3 Proof of Theorem [D 

We proceed to prove Theorem[T] Let / G H^HO, T) x dM) and u solve dOt . Put 
li' = dtu, then 

' (i(9t + Ag + g(x))u'(i,x) = in (0,r)xX, 

m'(0,x) =0 in (3.47) 
^ u'{t,x) = f{t,x) on (0,r)x5A^, 
Since /' G L2((0, T) x dM), by lemma[34l we get 

u' G C(0,T;i/-i(X))ni7-^(0,T;L2(7W)). (3.48) 
Furthermore there is a constant C > such that 

ll"'llc(0,T;//-i(A1)) + ll^1l//-i(0,T;L2(A4)) - ^ \\f\\m{{0,T)xdM) ' (^.49) 



Thus (13.48! ) implies the following regularity for u 

u G C^{{),T;H-^{M)) r\C{f),T;L^{M)), 
AgU G C(0, T; H-^{M)) n ^-^(O, T; ^^(A^)). 
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Since /(t, •) G H^{dA4), by the elliptic regularity, we get 

u G C(0,T; H\M)) n C\0, T; H-\M)). 
Moreover there exists C > such that the following estimates hold true 

ll^llci(0,T;/i'-i(A4)) - ^ \\f\\m{{0,T)xdM) ^ 



ll^g'"llc(0,T;//-i(A4)) - ^ \\f\\m{{0,T)xdM) ■ (3.50) 
Using Lemma [331 we find 

ll'"llci(0,T;//-i{7M)) - ^ \\f\\Hmo,T)xdM) ' 

ll'"llc(0,T;//i(7M)) - ||/||//i((o,T)xaA4) ' (3.51) 

The proof of dl.SI ) is as in Lemma [331 If one multiplies ( 11.21 ) by {N, Vgu)^, the 
arguments leading to (13.201 ) give now 

= [ [ idtu{N,Vgu)dv^dt+ [ [ Agn(iV,Vgil) dv"dt 
Jo Jm ^ JoJm 

^ LIm ^''"^^ "^""^ dt = I[+I'^ + (3.52) 

with 

|Re/{| < ||/||^i((o^y)^g_y^^) + e ||9jyn||^2((o,T)x9A^) ' (3.53) 

where we have used (13.511 ) instead of (I3.13I )- (I3.7I) . Furthermore, we derive from 
Green's formula 

Rel2 = l [ [ |a^n|^ dcj""Mt+ /" [ Z)giV(Vgn, Vgu) dv^ dt 
2 Jo JdM Jo Jm 

-If / |Vg7xPdiv(iV)dv^dt-i / / |V,/|2da^-Mt. (3.54) 
Jo Jm ^ Jo JdM 



This together with 

\m{{0,T)xdM) 



Re/^l < ||/|lHiao,T)xaA4) (3-55) 



and ([3351 ). (13341 ) and (l333l ) imply 

II^i'^IIl2{(o,t)xOM) ^ ll/ll/fi{(o,T)x9;n) ' (3.56) 
where we have used (13.51! ) again. The proof of Theorem [T] is now complete. 
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4 Geometrical optics solutions of the Schrodinger equa- 
tion 



We now proceed to the construction of geometrical optics solutions to the Schrodinger 
equation. We extend the manifold {Ai,g) into a simple manifold Ai2 3) Ai and 
consider a simple manifold (A^i,g) such that A^2 -A^i- The potentials (71,(72 
may also be extended to 7W2 and their H^{Aii) norms may be bounded by Mq. 
Since qi and q2 coincide on the boundary, their extension outside A4 can be taken 
the same so that qi = q2 in 7W2 \ -M-i. 

Our construction here is a modification of a similar result in ifTTI . which dealt 
with the situation of the wave equation. 

We suppose, for a moment, that is able to find a function G C^(A^) which 
satisfies the eikonal equation 

|V,*i:=E^«S|^ = l. V«A<. (41, 

l,j = l * ^ 

and assume that there exist a function a G H^{W, H'^{M.)) which solves the trans- 
port equation 

|+E8''||fe + 5(^.«-«. A. (4.2, 

j,k=l 

with initial or final data 

a(t,rc) = 0, VxGTW, and t < 0, or t > Tq. (4.3) 
We also introduce the norm ||-||^ given by 

IkL = l|a|l//i(o,To;//2(»)) ■ (4-4) 

Lemma 4.1 Let q G L°°[Ai). Then the following Schrodinger equation 

(i(9t + Ag + g(x))n = 0, in Mr ■■= {0,T) x M, 
u{k,x) =0, K = 0, or T 

has a solution of the form 

u{t, x) = a{2Xt, x)e^^(^(^)-^*) + vx(,t, x), (4.5) 
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such that 

u G C^{Q,T;L^{M))r\C{^,T;H'^{M)) (4.6) 
where v\{t, x) satisfies 

vx{t,x) = Q, \/{t,x) e {0,T) X dM, 
vx{k,x) = 0, X £ Ai, K = OorT. 

Furthermore, there exist C > such that, for any A > the following estimates 
hold true. 

\\Mtr)\\H^^M)<C\^~^\\a\l, /c = 0,l. (4.7) 
The constant C depends only on T and M ( that is C does not depend on a and X). 
Proof . Let us consider 



,x e 



iX{ip-Xt) 



(4.8) 



k{t, x) = - {idt + Ag + (a(2Xt, : 
Let V solve the following homogenous boundary value problem 

' {idt + /^g + q)v{t,x) = k{t,x) in (0,r)xA^, 

v{k,x) = 0, in M, T = 0, or T 

v{t,x) =0 on (0,r) X dM, 



To prove our Lemma it would be enough to show that v satisfies the estimates (14.71 ). 
We shall prove the estimate for k = 0, and the k = T case may be handled in a 
similar fashion. By a simple computation, we have 



(4.9) 



-k{t, x) = e^^(^(^)-^*) (Ag + q{x)) {a{2Xt, x)) 



j,k=i 



dip da a 



dxj dxj. 2 



+ -AgV' {2Xt,x) 



+ A^a(2At,x)e 



iX{4){x)-Xt) 



j,k=i 



dxj dxk 



(4.10) 



Taking into account (14.11 ) and (14.21) . the right-hand side of (14.101 ) becomes 

k{t, x) = _e*^('^(^)-^*) (Ag + q) {a{2Xt, x)) 
_ -e'^W^^-^i)kQ{2Xt,x). 



(4.11) 
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Since ko G H^{0,T; L^{M)), by Lemma[l2 we find 

vx G C\0,T;L^{M))nC{0,T;H\M)nH^{M)). (4.12) 
Furthermore, there is a constant C > 0, such that 

\\vx{tr)\\LHM)<C \\ko{2Xt,-)\\mM) (4.13) 

c r 

- ^ II II* 

Moreover, we have 

l|Vt;A(t,-)llL2(A1) 

<Cr/^ (A2||A;o(2At,-)|li2(^) + A||atA;o(2At,-)|li2(^)) dt 

+ ||A:o(2At,-)IL2(^) dt. (4.14) 

Finally, choosing r/ = A~^, we obtain 

l|VvA(i,-)llL2(Ai) < C' II^o(s,-)|Il2(A4) + jj\dtko{s,-)\\L^M) d's) 

<C||aL. (4.15) 

Combining (14.151 ) and ( 14.131 ). we immediately deduce the estimate ( 14.71 ). □ 

We will now construct the phase function ip solution to the eikonal equation 
(14.11 ) and the amplitude a solution to the transport equation (14.21 ). 



Let y € dAdi. Denote points in TWi by (r, 9) where (r, 9) are polar normal 
coordinates in A^i with center y. That is a; = e'x.py{r9) where r > and 

e G SyMl = G TyMl, l^lg = l} . 

In these coordinates (which depend on the choice of y) the metric takes the form 

g{r,9) = dr^ + go{r,9) 
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where go{r, 6) is a smooth positive definite metric. For any function u compactly 
supported in A^, we set for r > and 9 G SyMi 

u{r,9) = u{expy{r9)) 

where we have extended u by outside A4. An explicit solution to the eikonal 
equation (14.11 ) is the geodesic distance function to y G dA4 1 

V'(x) = dg(x,y). (4.16) 

By the simphcity assumption, since y G M-2\M., we have -0 G C°°{M.) and 

^[r,9) = r = dg{x,y). (4.17) 

The next step is to solve the transport equation (14.21) . Recall that if /(r) is any 
function of the geodesic distance r, then 

Ag/(r) = /"(r) + — — /'(r). (4.18) 

Here a = a{r, 9) denotes the square of the volume element in geodesic polar 
coordinates. The transport equation (14.21) becomes 



Thus a satisfy 



da ^ d'tjj da ^ ^aa~^ dadil: _ ^ 
dt dr dr 4 dr dr 

da ^da _ida ^ (4 20) 

dt dr 4 dr 
Let (j) G Cg°(M) and b G H'^{d+SM). Let us write a in the form 

a{t, r, 9) = a~^/V(t - r)b{y, 9). (4.21) 

Direct computations yields 

^(t,r,e) = a-i/V(t-r)%,e). (4.22) 

and 

^{t, r, 9) = -la-5/4^(/>(t - r)h{y, 9) - a-^l^^{t - r)h{y, 9). (4.23) 

Finally, (1423] ) and (l422l ) yield 

dH , ^, (9a, ^, 1 ^^da , . ^ 

— {t, r, 9) + — (t, r, 0) = --a-ia(t, ^)^. (4.24) 

Now if we assume that supp(0) C (0,eo)» > small, then for any x = 
ex.py{r9) G M, it is easy to see that a{t,r,9) = if t < and t > Tq for 
some To > sufficiently large. 
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5 Stable determination of the electric potential 



In this section, we complete the proof of Theorem |2] We aie going to use the 
geometrical optics solutions constructed in the previous section; this will provide 
information on the geodesic ray transform of the difference of electric potentials. 

5.1 Preliminary estimates 

The main purpose of this section is to present a preliminary estimate, which relates 
the difference of the potentials to the Dirichlet-to-Neumann map. As before, we let 
qi, q2 ^ =S(Mo) such that qi = q2 on the boundary dA4. We set 

q{x) = {qi - q2){x). 

Recall that we have extended qi,q2 as H^{M-2) in such a way that g = on 
M2\M. 

Lemma 5.1 There exists C > such that for any ai, a2 G //^(M, H'^{M.)) satis- 
fying the transport equation d4.2D with d4. JD . the following estimate holds true: 

nq{x)ai{2\t,x)a2{2Xt,x) dv" dt 
- .A 

<C(A-2 + ||Ag,,, -Ag,,,||)||ai||J|a2L (5-1) 

for any sufficiently large A > 0. 

Proof . First, if a2 satisfies (14.21) . (14.31) and A is large enough, Lemma l4T] guaran- 
tees the existence of a geometrical optics solution U2 

U2{t, x) = a2{2Xt, x)e^^(^(^)-^*) + V2,xit, x), (5.2) 

to the Schrodinger equation corresponding to the electric potential q2, 

{idt + Ag + q2{x))u{t,x) = m{0,T)xM, u(0,-)=0 inA^ 

where V2,\ satisfies 

A \\v2,xit, ■)\\l2^m) + l|V^2,A(t, < C \\a2\l (5.3) 

V2,x{t, x) = 0, \/{t, x) G (0, T) X dM. 

Moreover 

U2€C\0,T;L\M))nC{0,T;H^{M)). 
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Let us denote by fx the function 

fx{t, x) = a2(2At, x)e^^(^(^)-^*) , t e iO,T), x G dM. 

Let us consider v the solution of the following non-homogenous boundary value 
problem 

' {idt + Ag + qi)v = 0, {t,x) £ {0,T) X M, 

< v{0,x)=0, xeM, (5.4) 

^ v{t, x) = U2{t, x) := fx{t, x), (t, x) G (0, T) x OM. 

Denote w = v — U2- Therefore, w solves the following homogenous boundary 
value problem for the Schrodinger equation 

' {idt + Ag + qi{x)) w{t, x) = q{x)u2{t, x) {t, x) G (0, T) x M, 
< w{0,x) =0, X e M, 

w{t,x) = 0, {t,x) G (0,r) X dM. 

Using the fact that q{x)u2 G W^^^{0, T; L'^(M)) with u{0, •) = 0, by Lemma[321 
we deduce that 

w G C\0,T;L'^{M))nC{0,T;H^{M)nH^{M)). 

Therefore, we have constructed a special solution 

ui e C\0,T;L'^{M)) nC{0,T;H^{M)) 

to the backward Schrodinger equation 

iidt + Ag + qi{x))ui{t,x) = 0, {t,x) G (0,r) xM, 
ui{T,x) = 0, X e M, 

having the special form 

ui (t, x) = ai (2At , x)e^^(^(^)-^*) + vi^x x) , (5.5) 
which corresponds to the electric potential qi, where vi^x satisfies 

A \\vi,xit, ■)\\l2^m) + l|V^^i,A(i, -^mM) ^ C \\ai\l . (5.6) 
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Integrating by parts and using Green's formula (12.41 ). we find 

n{idt + Ag + qi) wui dv^ dt = [ [ g'U2«idVgdt 
„ .A Jo JM 



rT 
lO JdM 



duwuidag'^dt. (5.7) 



Taking ( 15.71 ). ( 15.41 ) into account, we deduce 

nq{x)u2{t, x)ui(t, x) dv" dt 
rT 



[ [ {Ag,,,-Ag,,,){fx){t,x)g^{t,x)da2-'dt (5.8) 

Jo JdM 



10 JdM 
wliere gx is given by 

gx{t,x) = ai(2At,x)e^^('^(^)-^*), (t,x) G (0, T) x OM. 

It follows from (|5^ . (1531 ) and (|5^ that 



/ / q{x){a2ai){2Xt,x)dv2dt = - [ [ gx{As,g, - Ag,g,) fxda^-' dt 
Jo Jm Jo JdM 

rp 

[ [ qe'^^^-^''^a2i2Xt,x)vi,xdv'2dt 
Jo Jm 



T 

qv2,xe-'^'-^~^'hii2Xt, x) dv" dt 

10 JM 

-II qv2.xvi,xd^ldt. (5.9) 
Jo Jm 

In view of (15.61 ) and (15.31 ). we have 

n ge*^W'-^*)a2(2Ai, x)v^,x dv" dt 
'0 Jm 

<C ll«2(2At,-)|li2(^)||t;i,A(t,-)IL2(^) dt 
<CX-^\\a2\lhi\l- (5-10) 

Similarly, we deduce 

rp 

I I ge-*^W'-^*)ai(2At,x)z;2,A(t,a;)dv"dt 
10 Jm 



<CX-^\\ai\l\\a2\l. (5.11) 
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Moreover we have 



{x)v2,\{t, x)vi^x{t, x) dv" dt 



10 JM 

On the other hand, by the trace theorem, we find 



< C\ ^ ||ai||^ ||a2||^ 



(5.12) 



JO JdM 



Ag,g,)(/A)5AdarMt 



< l|Ag,gi Ag,g2ll II/a||_h-1((0,T)xOM) II^A |Il2((o,T)x9A4) 



< c 



AV2 
Ai/2 



g, 91 Ag, 92 I 



(5.13) 



The estimate ([SlB follows easily from (ESI), (ISTTOl ). (ISHl ). (ISHl ) and (I57T3] ). This 
completes the proof of the Lemma. □ 



Lemma 5.2 Let Mq > 0. There exists C > such that for any b € H'^{d+SMi), 
the following estimate 




r+(yfi) 

q{s,e)b{y,e)ii{y,e)dsdu:y{d) 



<CU^,,,-!^^,,,f'\\b{v,-)U,^stM,)- (5-14) 

holds for any y £ dAii. 
We use the notation 

S+Mi = {eeSyMi:{u,e)g<0}. 

Proof . Following (14.211 ). we take two solutions to (14.21 ) and (14.31) of the form 

ai(t.r,e) = a'^/U{t-r)biy,9), 
a2{t,r,e)l = a'^^Uit - r)fiiy,9). 

We recall that fi{y,6) = {v{y),9) is the density of the space where the image of 
the geodesic ray transform lies. Now we change variable in (15.11 ). x = exp (r0). 
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r > and 9 G SyMi, we have 

nqai{2Xt, x)a2(2At, x) dv" dt 



T 




'0 JSyMiJO 
l-T 



q{r, e)ai {2Xt, r, 9)a2{2Xt, r, e)a^''^ dr dujy{e) dt 




JSyMiJO 
2XT 



1 

2A 



q{r, 9)4)'^ {2\t - r)b{y, 0)ii{y, 6) dr dujy{e) dt 

r+{y,0) 

q{r, e)4>^{t - r)b{y, 9)fi{y, 9) dr duy{9) dt. 



10 JSyMiJO 

By virtue of Lemma ISTTl we conclude that 



JSyMiJO 



r+{y,e) 



qir, 9)4>\t - r)b{y, 9)fiiy, 9) dr dojyi9) dt 



< C (A-i + A ||Ag,,, - Ag,,,||) ||0||^3(K) 11%, ■)\\mistMi) ■ (5-15) 

By the support properties of the function cf), we get that the left-hand side term in 
the previous inequality reads 



JSyMiJO 



r+(y,e) 



q{r, 9)(t>''{t - r)b{y, 9)fi{y, 9) dr dujy{9) dt 



ct>{t)dt 



{r,9)biy,9)f,iy,9)drdu;y{9). 



ISyMlJO 

Finally, minimizing in A in the right hand-side of (15.151 ) we obtain 



JSyMiJO 



-iy,9) 



{s,9)biy,9)f,{y,9)dsdu;yi9) 



< r' II A — A 11^/^ 

^ IM*-g, gi ^^g, 92 II 



\\Hy,-)\\H^(^s+Mi 



) • 



This completes the proof of the lemma. 



□ 
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5.2 End of the proof of the stabiUty estimate 

Let us now complete the proof of the stability estimate in Theorem |2l Using 
Lemma 15^ for any y £ dMi and b £ H'^{d+SM) we have 

/ Z{q){y,e)b{y,e)ii{y,e)du:y{d) 

JSyMl 

< C'll^g.qi -^g,q2ll^^^ll%.-)ll//2(5+^^). 

Integrating with respect to y G dM i we obtain 

/ Z{q){y,e)b{y,e){e,v{y)) dal^-\y,e) 
Jd+SMi 

(5.16) 

Now we choose 

b{y,e)=I{ri{q)){y,e). 
Taking into account ( I2.14I ) and ( I2.10I ). we obtain 

112^*2^(9) IIl2(A4i) ^ ^ ll^g,gi ~ ^g,g2ll^'^^ • 
By interpolation, it follows that 

\\^*AQ)\\m(M^) <C\\rimLHM^) \\^*^i<l)\\HHM^) 

< C ||X*X(g)||^2(_A4^) Ikll^HA^) 
<C\\I*imLHM^) 

<C||Ag,,, -Ag,,,f/^ (5.17) 

Using (12.131 ). we deduce that 

^ ^ ll^g,9i ~ ^g,g2ll^^^ • 
This completes the proof of Theorem |2l 

6 Stable determination of the conformal factor 

This section is devoted to the proof of the stability estimate for the conformal factor. 
We use the following notations; let c G ^{Mq, k, e), we denote 

go{x) = l-c{x), ^,i(x) = c"/2(x)-l, £>2(x) = c"/2-i(x)-l, 

q{x) = Q2ix) - Qlix) = c"/2-i(a;)(i _ c(x)). (6.1) 
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Then the following holds 



I|£'o|Il2(_a/,) < ||£'|Il2(_a/,) < C \\qq\\^2 



(M) 



(6.2) 



As in the case of potentials, we extend the manifold (7W , g) into a simple manifold 
M2 3) so that M-2 ^ Ml ^ M with (A^i, g) simple. We extend the confor- 
mal factor c by 1 outside the manifold 7W ; its C'^ {M. 1 ) norms may also be bounded 
by Mq. The first step in our analysis is the following lemma. 

Lemma 6.1 Let c £ C°°{Ai) such that c = 1 on the boundary dAi. Let ui, U2 
solve the following boundary problems in (0, T) x 7W with some T > 



^2(0,-) = 

U2 = /2, 

Then the following identity 



' {idt + Ag)iii = 0, in (0, T) x M 

ui(0, •) = 0, in M 

Mi = /i, o«(0,r)x(97W 
{idt + Acg)'U2 = 0, in (0, T) x M 

in M 

on (0, T) X dM 



(6.3) 



(6.4) 



JdM 



(Ag-Aeg)/i/2 da^-Mt 



Qi{x)uidtU2 dvgdt 



JM 



+ 



JM 



Q2{x) {VgUi{t,x),VgU2{t,x)) dvgdt (6.5) 



holds true for any fj G H^{{0, T) x dM), j = 1, 2. 

Proof . We multiply both hand sides of the first equation (16.31) by U2, integrate by 
parts in time and use Green's formula (12.41) to get 



= [ [ (iatui + AgUi)ll2dv"dt 
Jo JM 

rT 



T 



uidtU2 dv"gdt 



JM 



QiUidtU2 dv" dt 



JM 



and after using a second time Green's formula, we end up with 

Jo Jm Jo JdM 



j-T 

+ 

J{ 

+ f [ d,uj~2 da^-' dt. 

Jo JdM 



Taking into account tlie fact that c = 1 on dM, the fact that {—idtU2 + ^c£^^2) = 
in (0,T) X A1, and the fact that the Dirichlet-to-Neumann map is selfadjoint, it 
follows that 

ri (Ag - Aeg) /l ^ d(7^-l dt = i n Ql{x)uidtU^ dv^ dt 

Jo JdM Jo Jm 

This completes the proof of the Lemma. □ 



6.1 Modified geometrical optics solutions 

Let ^i, %l)2 be two phase functions solving the eikonal equation with respect to the 
metrics g and eg. 

onM. (6.7) 

IV7 / |2 jkd^2dll)2 -. 

3,k=l ^ ^ 
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Let a2 solve the transport equation in M x 7W with respect the metric g (as given 
in section 4) 

j,k=l ■' 

Let 03 solve the following transport equation in M x with respect to the metric eg 



j,k=i ^ 

= a2{t,x)ipoix,X) (6.9) 

which satisfies the bound 

llasL < CX WqoWc^m) II^sL • (6-10) 

Let us now explain how to construct a solution 03 satisfying (16.91) and (16.101 ). To 
solve the transport equation ( 16.91 ) and (16.101 ) it is enough to take, in the geodesic 
polar coordinates (r, 6) (with respect to the metric eg) 

a3it,r,e;X)=a-^/\r,9) [ all,\s,9)a2is-r + t,s,e)^o{s,9,X)ds, (6.11) 



where acg{r, 9) denotes the square of the volume element in geodesic polar coordi- 
nates with respect to the metric eg. Using that ||y'o("i '^)llci(A4) — ll^ollci(A^) 
and (16.111) we obtain (16.101 ). 

Lemma 6.2 Let c £ "^(Mq, k, e) be such that e = 1 near the boundary dM.. Then 
the equation 

(i^f + Acg) li = 0, in (0,r)xX, u(0,x) = (6.12) 
has a solution of the form 

U2{t,x) = ia2(2At,j;)e^^('''i(^)-^*) +a3(2At,2;;A)e^^('^2(x)-Ai) 

(6.13) 

which satisfies 

-^ll^^2,A(t,-)llL2(_^) + ||Vu2,A(i,-)llL2(^) + A"^ \\dtV2,x{t,-)\\^2(^M) 

<c(a2||£,o|Ici(m) + A-') \\a2L (6-14) 

where the constant C depends only on T and M ( that is C does not depend on a, 
X and e). 
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Proof . We set 

kit, x) = - (idt + Aeg) Qa(2At, x)e*^('^i-^*) - a3(2At, x, A)e^^(^2-Ai) 

To prove our Lemma it is enough to show that if v solves 

{idt + Acg)v = k (6.15) 
with initial and boundary conditions 

v{0,x) = 0, in M, and v{t,x) = on {0,T) x dM (6.16) 
then the estimate (16.141 ) holds. But we have 

-k{t, x) = ie*^^'^^-^*) Acga2(2At, x) (6.17) 
A 

+ 2^e-(^--) E cg^-^||g + f A.g^,j (2At,x) 

jf, din din 



+ XeiH^^-^^)a2{2Xt,x) 1 - ^ g- 
^giA(^.-Ai) (a3(2At,x)) 



dxk 



n 



+ A2e^A(^.-A*)„3(2At,x) I 1- Y.(^^y'^ — 

j,k=l ^ 



dxj dxk 



Taking into account (16.71 ) and (16.81) . the right-hand side of (16.171) becomes 
1 
A 



k{t, x) = -e*^^'^!-^*) Acga2(2At, x) (6.18) 



+2ie'H^^''t) l^^^-i _ 1) (VgVi, Vga2(2At, x))^ 
+ ^a2(2At,a;)(AcgV'i - AgV'i) 

j,A:=l 

+ ^e'^^'^^-^-^\l-c-')]{2Xt,x) 
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By ( |6!9l ) we get 

-k{t,x) = ^e^^(^i-^*)Acga2(2At,x) 
A 



+2ie*Al^i~-')(^(c-^-l)(VgV'i,Vga2(2At,x))^ 

+ ^a2(2At,x) (AegV'i - Ag-^i^ 
+ e^^(^2-^*)Aega3(2At,x) 

= ^e^^(^i-^*)A;o(2At,x) + e*^('^i-^*)A;i(2At,x) +e^^(^2"At)^^(2At,x). 
A 

Since kj e ifgHO^ L^{M)), by Lemma[l2l we deduce that 

vx €C\0,T;L^{M))nC{0,T;H^{M)nH^{M)) (6.19) 

and 



\\v\{t,-)\\L2(M) 







< C / - ||A;o(2At, + ||fci(2At, + ||fc2(2At, 



A 



dt 



<- I ( -||A:o(s,-)|Il2(^) + ||A;i(s,-)IL2(^) + ||/c2(s,-)IIl2(^) 1 ds 



A ./to V A 

< 



(7/1 

( ;^ l|a2L + II£'o|Ici(m) lk2L + A^ lle-ollcHM) ll^al 



< <^ (a||^o|Ici{m) + ^) ll«2L. 
Moreover, we have 

r-T 



\\^vx{t,-)\\LHM)<Cv\^ (A||A:o(2At,-)IL2(^) + ||at/co(2At,-)IL2(A4)) 

+ r (A2||A:i(2At,-)|lL2(^) + A||aifci(2At,-)|lL2(^)) dt 



+ r(A2||A:2(2At,-)|lL2(Ai) + A||atA;2(2At,-)IL2(^)) dt 







Ql|fco(2At,-)IL2(A4) + ll^i(2At,-)|li2(^) 

+ ||/c2(2At,-)|li2(^) Idt. 
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Hence, we obtain the following estimate 

l|Vt;A(t,-)llL2(A4) (6-20) 

+ ^ (a||A;i(s,-)|Il2(^) + \\0tki{s,-)\\L2(^M)) 

+ J (^xiMs 

+ ^ II^2(s,-)|Il2(a^)) dt 



')\\l^{M) 



<Cv(l + ll^ollci(Af)) l|a2L + " + ^ ll^ollci(M) 



02 



|a2|L. (6.21) 



Now choosing rj = \ ^, we obtain 

\NMt, ■)\\mM) <C(^J + X^ llf?o|lci(A/) 
Finally, if we study the equation satisfied by dtv, we also find 

\\dtv{t,-)\\L2iM)<c{l + X'\\Qo\\c^M)) \W2L. (6.22) 
This ends the proof of Lemma |6]2] □ 



Lemma 6.3 There exists C > such that for any ai, 02 G H'^{M)) satis- 

fying the transport equation (16.81 ) with (14. 3I ). the following estimate holds true 



A 



/" /" ^(a;)(aia2)(2At,a;)dv"dt 
'0 



< ll£'o|lci(A4) (a ^ + ll£'o|lci(Ai)j + A ||Ag - Acgll | ||aiL ll«2|L 

(6.23) 

for any sufficiently large A. 
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Proof . Following Lemma [6i2] let U2 be a solution to the problem {idt + Acg)u = 
of the form 

U2(t,x) = ia2(2At,x)e-*^('^i-^*) +a3(2At,x;A)e-*^('^2-^*) +V2,x{t,x) 
A 

where V2^x satisfies (16.141) and 03 satisfies (16.101) . Thanks to Lemma l4~T] let ui be a 
solution of {idt + ^g)^ = of the form 

ui{t,x) = ai{2Xt,x)e'^^^'~^^^ + vi,x{t,x). 

where vi x satisfies (14.41 ). Then we have 

dtMi,^) = 2dta2{'^Xt,x)e"^^^''~^'^ + iXa2{2Xt,x)e-'^'^^'~^'^ 

+ 2A9ta3(2At, x; X)e-^H-4'2-xt) ^ i^2-^(2At, x, A)e-*^('^2-At) 
+ dtV2,xit,x). (6.24) 



Let us compute the first term in the right hand side of (16.51 ). We have 
/ / QiUidtU2dv2dt = iX [ [ ^i(aia2)(2At,x)dv^dt + Ji(A) + J2(A) 

Jo JM Jo Jm 

(6.25) 

with 

Ji(A) = 2/ / ^i(ai5ta2)(2At,x)dv"di 
Jo Jm 

+ 2// QiVi^xdta2{'^Xt,x)e-'^^'^^-^'^dwldt 
Jo Jm 

+ iX [ [ a?;i,Aa2(2At,2;)e-^^(^i-^*)dv"dt 
JO Jm 



and with 



J2(A) = +2A /7 Si Mas) (2At, x)e'^^^'-^'^ dv" dt 
Jo Jm 

+ iX^ [ [ ^1(0103) (2At,x)e^^('^i-'^2)dv"dt 
Jo Jm 

+ Qiai{2Xt,x)dtV2,\ii,x)dVgdt 
Jo Jm 

+ 2X f f QiVi^xdtM'^Xt, x)e-*^(^2-^*) dv" dt 
Jo Jm 



38 



JM 



1 

10 JM 

^From (16.101 ). ( I6.14I ) and ( |4.4b we have the estimates 



nQiVi,\dtV2,\{t,x)d\'i dt. 
_ . A 



\JiW\<C\\q4cHm) A"'ll«2||J|ai|L 

|J2(A)| < C WqoWc^m) (^"' + ll^?o|lci(A4)) ll«2L ll«iL • (6-26) 
On the other hand, we have 

Vgui = (Vgai)(2At,x)e*^('^i-^*) + iA(Vg^i)ai(2At,x)e*^('^i-^*) + VgZ;i,A 

Vgll2 = Y(Vga2)(2At,2;)e-^^('^i-^*) - m2(2At, 2;)VgVie-^^('^i-^*) 
A 

- iAa3(2At,x)VgV'2e"^^('^^"^*) + V ^a?X2\t,x)e-'^^'^-'-^^^ + VgU2,A 
and the second term in the right-hand side of ( I6.5I ) becomes 

^ ^ ^2(2;) (Vg'Ui(t, x), VgU2(t, x))g dv^ dt 

= x[ [ e2(x)(aia2)(2Ai,x)dv"dt + J3(A) + J4(A) (6.27) 

Jo JM 



with 

J3(A) 



t// ^>2(x)(Vgai(2At,x),Vga2(2At,x)) dv"dt 
X Jo JM 
T f 

Q2ix)a2{2Xt,x) (Vgai(2At,x),VgV'i(x))^ dv" df 
10 JM ^ 

+ i Q2{x)ai{2Xt,x){V„a2{2\t,x),VMx))dyldt 
Jo JM 

+ \ I I £^2(x)e-*^(^i-^*)(Vga2(2At,x),Vg^;i,A(t,x)) dv^dt 
^ Jo JM 
T r 

Q2ix)a2{2Xt,x)e-'^'^^'-^'^ {V^v^,xit,x),V^Mx))^ dv" dt 
'0 JM ^ 

and with 

J^{X) = -i\ [ [ £>2(x)a3(2At,x)e*^('^i-'^2)(Vgai(2At,x),VgV'2(x))^ dv"dt 
Jo JM 
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+ /V ^?2(x)e^^(^i-'^2)^Vgai(2At,x),Vga3(2At,x)) dv"dt 
Jo Jm 

rT 



+ / f?2(x)e^-'^'^^-'^'^(Vgai(2At,x),Vg?J2,A(t,x)) dv^dt 
JO Jm 

LIm ^2(^)("i"3)(2At,x)e*^('^i-'^2) (Vg^i(x), VgV2(x))g dv-dt 

+ iX [ [ ^>2(x)ai(2At,x)e*^('^i-'^2)^y^^3(2At,x),VgV'i(2;)) dv"dt 
Jo Jm 

+ iX [ I ^2(x)ai(2At,x)e*^('^i-^*)(Vg^i(x),VgU2(t,x)) dv"dt 
Jo Jm 

-ill i?2a2(2At,x))e-^^(^i-^*)(Vgi;i,A(t,x),VgVi(x)) dv^dt 
- iA /" /" ^>2a3(2At, x)e-^^('^2-^*) (Vg?;i,A(t, x), VgV2)g dv^ dt 

T r 



+ / / f?2e-*^^'^^-^'^(Vg7;i,A,Vga3(2At,x)) dv^dt 
JO JA1 

+ / / f?2(Vgt;i,A,VgI;2,A>„ dv"dt. 
JO J7M 

^From (16.101 ). (16.141 ) and (l44l) . we have 

l^3(A)| < ||^>o|lci(M)A^M|a2LI|aiL (6-28) 

\JiW\ < ll^ollci(AI) (^"^ + A^ ll^ollci(A4)) II«2|L lloilL • 

Taking into account (I6.5I ). (16.251 ) and (16.271) . we deduce that 

/V (^g - ^-g) /i72 df^r ' = A Z"^/ i5(x)(aia2)(2At, x) dv" dt 
Jo JdM Jo Jm 

+ MX) + J2W + J3(A) + Ja{\). (6.29) 
In view of (16.261 ) and (16.281 ). we obtain 



A 



Q{x){aia2){2Xt^ x) dvg dt 
'0 JM 



< 



C{ ll£'o|lci(M) (a ^ + II£'oIIci{M)) + A ||Ag - A^gH } ||ai||, ||a2| 



This completes the proof. □ 
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6.2 Stability estimate for the geodesic ray transform 

Lemma 6.4 There exists C > such that for any b G H'^{d+SM-i) the following 
estimate 



I{Q){y,9)b{y,e) {9,u{y)) daf-^(y,0) 



/ 

Jd+SMi 

< C(^{\ ^ + ||£'o|lci(M) ) ll£'o|lci{M) 



+ A llAg - A, 



eg I 



\H^{d+SMi) 

(6.30) 



holds for any y € dMi. 

Proof . Following (14.211 ). we take two solutions of the form 

ai{t,r,e) = a-^/Uit-r)b{y,e), 
a2{t,r,e) = a~^/U{t-r)fi{y,e). 

Now we change variable in ( 16.231 ). x = e'x.py{r6), r > and 9 G SyM-i. Then 



n0ai{2Xt, x)a2{2\t, x) dv" dt 
_ .A 



JSyMiJO 
T 



-(yfi) 



JSyMiJO 
2XT 



1 

2A 



JSyMiJO 



>{r, 9)ai (2At, r, e)a2 (2At, r, Q)a^l'^ dr dbjy(B) di 
Q{r, 9)4? (2\t - r)b{y, 9)fi{y, 9) dr dujy{9) dt 

Q{r, 9)(j)\t - r)b{y, 9)fiiy, 9) dr du:y{9) dt. 



We conclude that 



T 

JSyMiJO 



r+{yfi) 



Q{r, 9)(t>\t - r)b{y, 9)fi{y, 9) dr dw^(0) dt 



+ A||Ag,g, -Ag,,, II ) 11^^2(5+^,) (6.31) 
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where SyMi = {9 ^ SyM : {0,iy)g < 0}. Given the support properties of the 
function (p, the left-hand side of the inequality reads in fact 

Q{r, e)4>\t - r)b{y, e)fi{y, 9) dr dojy{e) dt 



oo J SyMlJO 



^'{t)dt\x / Q{r,e)b{y,e)fi{y,e)drdojy{9). 

Integrating with respect to y G dA4 1 in (16.311 ) we obtain 

Iig){y,e)biy,e){e,,.{y)) dal^-\y,e) 

d+SMi 

cgll ) \\b\\H2(d+SMi) ■ 

This completes the proof of the lemma. □ 



6.3 End of the proof of Theorem |3] 

Let us now prove Theorem [3] We choose 

b{y,e)=I{I*Iiq)){y,9) 
and obtain using Lemma |64] and (12.101 ) 

0o\\c^(M) ) ll£'o|lci(A4) 

By interpolation we have 

\\'^*AQ)\\miM,) < C|l^*^(^?)llL2(A4,) \\^*Aq)\\hhm,) 
<c({\-'^ + \^\\eo\\cHM))\\Q4c^(,M) 

+x\\Kg- k^S)"^ \\x*x{e)\\%\j^y 

We use (12.131 ) and (12.141) to deduce 

1 g 

ll^lli2(A4) ^ ((^"^ + ll^o|lci{M) ) ll^o|lci{M) + ^ - ^cgll ) ' 

< + A2 ||^>o|Ici{^) )^ \\Q4h{M) + ll^'ollli(^) ||Ag - Aegll^ 
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Minimizing A ^ + 




1/2 



-eg 



1/2 



Since 



£'ollci(A^) 



< 



£'0||_H-n/2 + l + e(_V() 

,,24/25 II ||l/25 



< 



< 



1 1 24/25 
Q0\\l2{M) 



we therefore obtain 



6\\l'2{m) 



1/2 



-eg 



Taking e > small enough we conclude and obtain ( I1.12I) . 
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